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Abstract With the reliance in the processing of quantum 
information on a cold trapped ion, we analyze the entangle- 
ment entropy in the ion-field interaction with pair cat states. 
We investigate a long-living entanglement allowing the in- 
stantaneous position of the center-of-mass motion of the ion 



to be explicitly time dependent. An analytic solution for the 
system operators is obtained. We show that different nonclas- 



sical effects arise in the dynamics of the population inversion, 
■ depending on the initial states of the vibrational motion. We 
study in detail the entanglement degree and demonstrate how 
the input pair cat state is required for initiating the long liv- 
ing entanglement. This long living entanglement is damp out 
with an increase in the number difference q. Owing to the 
properties of entanglement measures, the results are checked 
using another entanglement measure (high order linear en- 
tropy). 



1 Introduction 

Controlling quantum systems and the couplings between 
various quantum degrees of freedom is one of the major 
goals in quantum information processing 1 . In the early 
days of quantum optics, the dynamics and behavior of 
typical quantum optical systems, such as the laser, was 
dominated by dissipative effects and uncontrolled fluctu- 
ations of system parameters. Systems of single trapped 
ions PlElEl and single atoms strongly interacting with 
a high-Q cavity were developed, corresponding to single 
quantum systems which can be prepared microscopically 
and observed under conditions close to idealized and fun- 
damental theoretical models |5j . By combining the tech- 
nologies of ion trapping and cavity quantum electrody- 
namics, deterministic coupling of single ions to an optical 
field has been achieved 

There is currently a diverse array of research interests 
being united under the banner of quantum information 
science, with fundamental research, in both experimental 



2 

and theoretical arenas, providing new insights on differ- 
ent time dependent interaction models [ZIIHIEI- Of great 
interest to physicists is the renewed vigor with which 
questions regarding the basic understanding and inter- 
pretation of quantum entanglement are being pursued. 
Also, the development of quantum information science 
and technology carries the potential of revolutionary im- 
pact on many aspects ^Hli as evidenced already by the 
applications in quantum cryptography, quantum com- 
munication, and rudimentary quantum computing. En- 
tanglement is today considered a fundamental resource 
in nature when it comes to quantum computation and 
information 11 , and measures of entanglement has be- 
come a major field of research \12\. 

Despite much effort and spectacular advances from 

long surviving entanglement remain elusive in the ion- 
field based systems. It is in this light that this paper 
will focus on long- living entanglement. The main moti- 
vation for the present work is its relevance to the field 
of quantum information, which is attracting broad in- 
terest in view of its fundamental nature and its poten- 
tially revolutionary applications. Our main purpose is 
twofold, first we wish to demonstrate how a pair cat 
state affect the entanglement for the ion-field interac- 
tion, and second we wish to see how the time-dependent 
amplitude of the irradiating laser field affect this entan- 
glement. The improved understanding of entanglement 
has attracted much attention recently from the quantum 



M. Abdel-Aty 

information community which may lead to new insights 
into the physics of correlated systems in quantum optics. 
Therefore there is a great need to compute these in an 
efficient way. In particular, we consider a single trapped 
ion which can be laser cooled to the ground state of 
the trapping potential and discuss the roles played by 
the initial state setting and time-dependent amplitude 
of the laser field on the entanglement. 

The paper is structured as follows: Section 2 is de- 
voted to a brief description of the pair cat state and 
show that the quadrature distribution reflects the en- 
tanglement. In section 3, we present the model of the 
completely quantized system and obtain an exact an- 
alytical solution of system operators. In section 4, we 
briefly discuss those features concerning the atomic in- 
version, which are relevant for collapse-revival phenom- 
ena. In section 5, we analyze in detail how the pair cat 
state affect the general features of von Neumann entropy 
(a measure of the ion- field entanglement). This gives us 
the opportunity to stress the essential role played by 
pair cat states in this context and study the existing of 
a long living entanglement. We summarize our results at 
the end of the paper and make some conclusions. 

2 Pair cat states 

Superpositions of orthogonal states which exhibits macro- 
scopic features appear to be of fundamental importance 
in recent studies of the foundations of quantum mechan- 
ics pni^ . Superposition refers to a quantum system 
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existing simultaneously in multiple states. This is usu- 
ally told as the parable of Schrodinger's cat, shut up in 
a box with a vial of cyanide that at any moment might 
be triggered to release its deadly gas by a radioactive 
decay reaction. Resonant microcavities can be used to 
study the behavior of mesoscopic superposition coher- 
ent states. Some success has been reported in creating 
such superposition states within high Q cavities in the 
optical domains \2'6\ . 

Recently, schemes have been presented for the prepa- 
ration of two-mode motional states of a trapped ion, 
such as pair coherent states ^j, pair cat states |25) . 
two-mode SU(1,1) intelligent states [201, and SU(2) cat 
states [22]. These schemes operate in two-dimensional 
isotropic traps. In terms of conventional coherent states 
parameterized in terms of a complex number a |2l| , the 
Schrodinger cat state has the following form 



yjl + 2cos(/)exp(-|ap) 



(la)-He^l-a)) 



(1) 



where ja) is a coherent state of amplitude a, and </) is a 
real local phase factor. Note that the relative phase can 
be approximately controlled by the displacement opera- 
tion for a given cat state with a >> 1 |28II29| . A special 
set of coherent states of the Barut-Girardello type known 
as pair coherent states has been formulated . If (a) 
and t'^(fe) denote two independent boson creation (an- 
nihilation) operators, then a6(a^5t) stands for the pair 
annihilation (creation) operator for the two modes. The 






Fig. 1 Contour plot of the quadrature distribution function 
P(x,y) for the pair cat state for g = 0, = 7r/2 and for 
different values of the parameter ^ where, (a) ^ — 0.1, (b)^ = 
1 and (c) 5 = 3. 
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pair coherent states \^,q) are defined as eigenstates of 
botli the pair annihilation operator ab and the number 
difference operator a^a— b^b , i.e., 

ab\tq)=m,q), [a^a~Vb)\tq)^q\^,q) (2) 

where ^ is a complex number and q is the charge pa- 
rameter, which is a fixed integer. Furthermore, the pair 
coherent states can be expanded as a superposition of 
the two-mode Fock states, 



\i,q)^N,Y, 



„=0 VnHn + qy. 



\n,n + q), (3) 



where Nq is the normalization constant and can be writ- 
ten as Nq — ''^(j(2|Cl)] ^ ; where Jq is the modified 
Bessel function of the first kind of order q. Pair coherent 
states are regarded as an important type of correlated 
two-mode state, which possess prominent nonclassical 
properties such as sub-Poissonian statistics, correlation 
in the number fluctuations, squeezing, and violations of 
Cauchy-Schwarz inequalities |31II32| . The pair cat states 
are defined as a superposition of two pair coherent states 
separated in phase by 180° 



N^{\^,q)+exp{t<j>)\~^,q)), (4) 



where the normalization factor is given by 



N^^ — (l + N^co8<pJ2 



n\(n + g)! 



(5) 



It has been shown that these cat states may exhibit 
stronger nonclassical effects than the corresponding pair 
coherent states. The correlated two-mode cat state |^, q, (j>) 
is the eigenstate of the operators a^6^ and a^a — b^b 






Fig. 2 Contour plot of the quadrature distribution function 
P{x,y) for the pair cat state for ^ = 1,0 = 7r/2 and for 
different values of the parameter q where, (a) q — 2, {h)q = 5 
and (c) q = 10. 
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with eigenvalues and q, respectively. In order to gen- tributions of the different components of the pair cat 

erate such cat states for the motion in a two-dimensional states are located close to the phase space origin, com- 

anisotropic trap we require two laser beams. Due to the peting with each other. When the values of the parame- 

strong nonclassical nature of the two-mode cat states, ter ^ is increased, the shape of the quadrature distribu- 

the generation of such states is of interest in testing tion becomes more pronounced by involving multi-peak 



structure (symmetric peaks). This is in a good agreement 
with the general behavior of the quadrature distribution 
of the pair coherent state |34| . It is interesting mention- 
ing here that, the two-mode pair cat states exhibit strong 
nonclassical features due to the correlation between the 
where {x\n) is a harmonic oscillator wave function. In two modes. 



quantum mechanics. 

The coordinate space wave function is given by 

{x,y\S) = N^{{x,y\£„q) + exp{i(l)){x,y\ -^,q)) 

OO 

= ^ (1 + (-l)"e»^) {x \n + q){y\n) , (6) 

11=0 



this case, the quadrature distribution is given by 

(1 + (-l)"e*'^) 2-(2«+9) 



Pix,y)^\N^\^ 



ri=0 



X Hn+qix)Hniy)e 



-0.5{x: +y ) 



n\(n + q)l 

2 



(7) 



In figure 2, we plot the quadratic distribution for dif- 
ferent values of the number of quant q, (say q — 2, 5 and 
10). We have shown that the behavior of the quadratic 
distribution makes possible, for sufhcient small values of 



For fixed values of the parameter q the periodic nature the indicator parameter q, to transform the multi-parts 
in the quadrature distribution function in the different windows to a one circular part. This is a consequence 
values of the parameter ^ are shown in figure 1, where, of the intrinsic properties of the pair cat state. Never- 
we plot P{x,y) for ^ = 0.1, 1 and 3. It is shown that theless, perfect purification is not achieved in this case, 
the regions of the quadrature distribution in the (x, y)- compared with the case in which small value of ^ is con- 
plane are symmetric with respect io y = —x and with sidered (see figure la and 2a). In figure la the window is 
respect to y = x (figures lb and Ic). It is clear that the not a two parts function, but is very narrow so that the 
present distribution refiects the entanglement present in conditional moment is approximately equal to one circle. 



the two-mode pair cat states. The main curiosity is that, 
if squeezing is strong enough, the ring- like quadratic dis- 



this behavior may be obtained if the parameter q takes 
smaller values. Note that the quadratic distribution is 



tribution function collapses to a quasi one dimensional dramatically different from the case when we consider 



object with a cigar form (see figure la). 



large values of q (see figure 2a and 2b), and indeed here 



The origin of such behavior is in the generation of it is not possible to observe separate parts of the distri- 



pair cat states in the microscopic regime, where the con- bution. It is interesting to show in this figure that the 
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Fig. 3 Contour plot of the quadrature distribution function 
P{x, y) for the pair cat state for g = 2, = and for different 
values of the parameter ^ where, (a) ^ = 1, (b)^ = 3 and (c) 



regions of the quadrature distribution in the (x, ?/)-plane 
are symmetric with respect to y = and with respect to 
X = (figures 2b and 2c). 

In figure 3, we set the parameter = and take 
different values of q and ^. In this figure we show the 
evolution of y) and its responses in a type of plot 
similar to that shown in figures 1 and 2. It is shown that 
the regions of the quadrature distribution are symmet- 
ric with respect to a; = and with respect to y = 
(figures 3b and 3c), but the general feature of this case 
is quit complicated compared with the previous cases. 
This can be understood from looking to equation 
where, when we set = this allows the contribution of 
the second term in equation jsj- Some of the differences 
are seen in the contour interval if a very small change of 
the controller parameters (^ and q) is considered. This 
can be demonstrated from a contour plot, such as shown 
in figure 3. 

3 Physical system and model Hamiltonian 

In the interest of retaining as much clarity as possible, 
we first recall some well-known facts about the system 
we wish to treat here, which consists of a single trapped 
ion interacting with a laser field 2T . We consider the sit- 
uation where oscillating classical electromagnetic fields 
propagating in two-dimensional trap irradiate a single 
trapped ion. To see that let us consider the laser fre- 
quency is tuned to the frequency difference between the 
two modes of the trap. Consequently the Hamiltonian 
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Using equation (10), the interaction part of equation (8) 
can be written as 



for such a system is given by P^25ll2()j . 

H — Hp + Ha + Hin , 



(8) 



where 



Hp = Tiujia\ai + Tiuj2a\a2, 
Ha = —^z, 



Hin = -p-E = -p.e{vt)[e^p[i{kix + k2V - ujt))] + H.c], 

with p the dipole matrix element of the ion p = p'^(ct+ + 
tT+), £{vt) is a modulated amphtude of the irradiating 
laser field. We denote by Wi, (i — 1,2) the fields frequen- 
cies, wq the natural frequency of the ion, and v denotes 
the ionic velocity. Here the a' s are the usual 2x2 Pauli 
matrices satisfying 



[o'2,o-±] = ±2(T±, [(T+,(J_] = 



(9) 



We denote by and a| the Bose operators for the quan- 
tized field mode which obey [a|, Oj] — Sij, where Sij — 1 
if i — j and otherwise. The operators x and y are the 
center of mass position of the ion which are quantized 
in the forms x — Ax{a\ + oi), y = Ay{a\ + 0,2) and 
Ax = v^l/(2wa;A/), Ay = ^y^/{2vyM) with M being 
the mass of the trapped ion. 

Making use of the a special form of Baker-Hausdorff 
theorem 35 the operator exp[i77(at + a)] may be written 
as a product of operators i.e. 



^iri{a'^ +a) 



exp 



?7^.t -1 



exp {irjo)^ exp [irjo) 



^e-^'/^^^M-^ ^ (!^)_^. (10) 



Hi; 



~p".e{vt)riir]2 exp 



1 



EE 

/=0 m=0 



(l + l)!(m + l)! 



(«'7i)'(*fy2)' 



(11) 



with 771 = kiAx and 772 = k2Ay. When we take Lamb- 
Dicke limit and apply the rotating wave approximation 
then the effective interaction Hamiltonian (10) takes the 
two photon bimodal interaction form. However, provided 
we adjust the strength of the electric field to include the 
time dependent factor e{vt). In this context we may refer 
to the work given in Ref. [22]. In this work the authors 
considered the usual JCM taking into account the ef- 
fects of stochastic phase fluctuations in the atom-field 
coupling coefficient. This is achieved by modulated the 
coupling parameter to include the time factor, in this 
case it is observed an interesting feature such as the de- 
coherence effect in the collapse and revival phenomenon 
of the Rabi oscillations. 

In the Lamb-Dicke regime, rjj << 1, equation Hll|l . 
can be well approximated by expanding the exponential 
terms up to second order in 77 (ji — Vi = ^2), 



Hin — hX{t)a\d2CT- + H.c, 



(12) 



ml 



where A(t) = — e(ut)p"77^ exp (^~^) ■ In what follows, 
we will find the solution of the equations of motion for 
the operators ai{t), (T±(t) and az{t) in the Heisenberg 
picture. Therefore the equations of motion related to the 
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Hamiltonian (8) can be written as 



4 Atomic inversion 



da I 
~dt 
da2 
~dt 
da- 
HT 
daz 
~dt 



-iuj2a2 — iA(t)ai(T-|_ 



(13) 



2i 



\{t) {d\a2cr 



- — aia2cr+ 



After straightforward calculations, one can find an ana- 
lytic time dependent solution as 

ai(i) = e*'^"(*) ( cos(7ia(t)) - ^ sin(7ia(t)) ) ai(0) 



-e-^'^^^e''''^^''^ sin(7ia(i))a2(0)<7_(0), (14) 



a2{t) 



71 

—iuj2ti{ia.{i) 



cos{j2Ct{t)) sin(72a(t)) ) a2(0) 

72 



72 



-»^.tg»A"(t) sin(72a(t))ai(0)CT+(0), (15) 



(16) 



(17) 



o-_(t) = _J_e-"^i2*e»'^"(*) sin(73a(t))ai(0)4( 

73 

+e-'-''''e'^'''^'^ (cos(73a(i)) 

+ ^sin(73a(t)) ) a-iO), 
73 / 

azit) = cos(274a(i))(Tz(0) + — sin(274a(i)) 

74 

X (^d\aia+{0) — a|a2(7_(0) 



where Wi2 — toi — lu2, 



73 = (^-^)(™+^), 74 = (^ + ^)(m+i), 

a{t) — / X{t )dt , /x = (a|a2<5'- + a2aiO'+). 
Jo 



Having obtained the dynamical operators, we are there- 
fore in a position to discuss some statistical properties 
of the system. 



The entanglement of motional degrees of freedom of the 
center of mass of the ion with it's internal degrees of free- 
dom manifests itself in the well known collapse and re- 
vival of population inversion ■ We shall start with the 
atomic inversion from which we can discuss the collapse 
and revival phenomenon. The (internal level) ionic dy- 
namics depend on the distributions of initial excitations 
of both the field and the center-of- mass vibrational mo- 
tion, given by {n\pf{Q)\n) = /3^„(0) and (m|p^,(0)|m) = 
Pmm(0)i respectively. For instance, the atomic popula- 
tion inversion may be written as 



(a,(t)) =Tr(a,(0)p(0)). 



(18) 



that is the difference between the probability of finding 
the system in the ground state and the probability of 
finding the system in the excited state. 

To analyze the effects resulting from variation in the 
parameter \[t) on the atomic inversion we have plotted 
in figure 4 the function {<Jz{t)) against the scaled time 
Xt, for different values A(f), where \{t) = A for figure 
(4a) and \{t) ~ Asinh(n7t) for figure (4b), where vj is 
a positive integer. In this figure let us take the param- 
eters f — 10,(7 = 1 S'lid (j) = 7r/2, keeping in mind we 
are dealing with the problem at exact resonances. In fig- 
ure (4a) we find that the value of the atomic inversion 
decreases from the maximum to its minimum, (approx- 
imately —0.75) then it collapses, but then it starts to 
fluctuate around zero for a short period of time. Thus 
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Fig. 4 The atomic inversion as a function of the scaled time 
At. Parameters: ^ — 10, g = 1, = 0, where (a) A(t) = At and 
(b) A(t) = Asinh(rot) 



the Rabi frequencies become commensurate so that the 
atomic inversion has an exact periodic evolution with 
a scaled A-dependent period. In this case there is no 
oscillations in the collapse region, as shown in figure 
4a. A different pattern for the oscillations occurs if the 
number difference is slightly changed. Strictly speaking, 
when we consider the time dependent modulating func- 
tion X{t) ~ Asinh(-n7t), the function sinh(n7t) oscillates 
(however their period of oscillation is not constant) then 
the appearance of the collapses and revivals in almost a 
periodic way has been washed out in this case. On the 
other hand and in the presence of the time-dependent 
modulating function we find the behavior of the atomic 
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inversion function is drastically changed where we can 
see irregular collapses and revivals occur, (see figure 4b). 
The short time revivals are strongly suppressed, as we 
see in figure 4b. We may compare the atomic inversion 
in both figure 3a and 3b, in the former case the collapse 
time is long enough to still allow revivals, while in the 
latter a shorter collapse time with the same amplitudes 
of the revivals. 

It is interesting to explore to what conditions can be 
considered to generate a long living entanglement. This 
will be seen in the next section. 

5 Entanglement 

It is well known that entanglement between several par- 
ticles is the most important feature of many quantum 
communication and computation protocols. In order to 
get a feel for how the entanglement of ion-field interac- 
tion affected by the presence of pair cat states we began 
by studying the aforementioned system that is analyti- 
cally solvable and start from a factored initial state of 
both the ion and field. We consider the notion of entan- 
glement entropy, that is, considering a quantum system 
in a pure state with the density matrix of the remaining 
space as pA — Trp [p). The density matrix contains all 
information of any system in a mixed or pure state, and 
computing vital physical information on any such system 
is mostly determined by the eigenvalues of the density 
matrix. The partial entropy of entanglement for a bipar- 
tite pure state is defined as Von Neumann entropy of the 
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reduced state, 

SA = -Tr{pA\og^PA}. (19) 

This procedure is well established, and works well for all 
cases where the initial state of the system is in a pure 
state. Most work has been focused on this entanglement 
in time-independent modulating function, but we will 
focus on a time-dependent case. 

On the other hand, the quantum entropy of the field 
Sf{t) can be expressed in terms of the eigenvalues X^{t) 
of the reduced field density operator as follows 

Sf (t) - - [A+ (t) In A+ (t) + XJ (t) In XJ (t)]. (20) 

In the case of a disentangled pure joint state Sf{t) is zero, 
and for maximally entangled states it gives In 2 [201 . Due 
to the higher dimensionality of the problem we cannot 
obtain a simple analytical expression for equation H20|) . 
therefore the numerical approach becomes indispensable. 

Now let us discuss the behavior of the field entropy 
resultant of the existence of the time dependent modu- 
lating function X{t). Figure 5 is a plot of quantum field 
entropy Sf{t) as a function of scaled time parameter Xt 
for 77 = 0.2, A (t) = A and considering different values of 
where ^ = 10 for figure 5a and ^ = 20 for figure 5b. Our 
numerical analysis was performed using other parame- 
ters from recent experiments |23II37) . where a trapped 
^Be"*" ion was laser-cooled to the zero-point energy and 
then prepared in a superposition of spatially separated 
coherent harmonic oscillator states. In |2H1 there have 
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Fig. 5 The quantum field entropy Sf{t) as a function of the 
scaled time Xt. The parameters are q = 1, r] — 0.2, \{t) = A 
and different values of ^ where (a) ^ = 10 and (b) ^ — 20. 

been recent demonstrations of the coupling between cav- 
ities and single-ion traps, where, a single calcium ion is 
stored in a spherical trap and placed in the center of a 
confocal resonator. 

It is interesting to note that, the maximum entangle- 
ment is decreases as the parameter q increased (see figure 
5a). As we see from figure 5a, as time goes on we note a 
growth in Sf{t), followed by a sudden decrease, almost 
down to zero at large values of q. When the time in- 
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creases further we see that the entanglenient has a long 
surviving, if we keep smaller values of q. The constat 
value of Sf{t) is interpreted as a result of quantum en- 
tanglement between ion and fields, surviving for a long 
time interaction. However, an interesting situation may 
arise when the parameter ^ takes large values (say 20 
in figure 5b). In this case, for some small fixed values 
of q, the entanglement surviving starts at earlier time. 
However, the overall features of the process are similar in 
both small and large values of ^. An important difference 
occurs only at earlier interaction time, where the num- 
ber of oscillations is increased with increasing ^. Also, 
the amplitudes of these oscillations are bigger compared 
to the case in which ^ = 10 (see figures 5a and 5b). As 
seen in the entropy plot in figure 5b, the pure state is 
basically obtained for large values of the parameter q 
(because the natural scale of q is between and ^, in fig- 
ure 5b, we expand the scale of q to take values between 
and 20). 

On the contrary, as can be seen from figure 5, if 
the fields are initially prepared in a pair cat state and 
the modulating function is taken to be time-dependent 
X{t) — Asinh(ti7t), then the entanglement parameter is 
significantly larger than zero for any time t > 0, but with 
some oscillations from time to time. Also, the number of 
oscillations is increased as the parameter q decreased as- 
sociated with low entanglement reaches zero for smaller 
values of the parameter ^ (see figure 6). 
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a 




Fig. 6 The same as figure 5 but A(t) — Asinh(cc7t). 

Another interesting aspect is that the time-dependent 
modulating function A(t) has an important affect on the 
entanglement, where, the total ion-field state can not 
have its purity diminished, which means that as the field 
becomes more pure the ionic state must be closer to a 
mixed state. Also, in this case, it seems remarkable that 
entanglement survive for long time but some oscillations 
occur here. More surprisingly still, survive is long-lived 
for the time-independent case. We should pointed out 
that if we simply use a time-independent modulating 
function, we would expect entanglement to survive. 



12 

Practically speaking, as can be deduced from figure 
5a, the oscillations in degree of entanglement between 
the ion and the field quickly damp out with an increase 
in q. The subsystem will not disentangle from each other 
until the steady state is reached as long as the interaction 
is present. This means that neither the ion nor the field 
will return back to a pure state except at the beginning 
of the interaction and t —^ oo, contrast to the q — case. 
In order to gain insight into the general behavior we have 
considered different values of ^ in figure 6b. The direct 
comparison of figure 6a and figure 6a, the amplitude of 
the oscillations in the degree of entanglement is increased 
with increasing the parameter f . 

The remaining task is to identify and compare the re- 
sults presented above for the entanglement degree with 
another accepted entanglement measure such as the lin- 
ear entropy. The question of the ordering of entangle- 
ment measures was raised in reference |39| . It was proved 
that all good asymptotic entanglement measures are ei- 
ther identical or fail to uniformly give consistent order- 
ings of density matrices jlH]. The best understood case, 
not surprisingly, is the simplest. In order to better char- 
acterize the subsystems, we may refer here to another 
measure of the entanglement of a reduced density ma- 
trix which is the product state identification 

Sl = 1- Trp% n > 2, 

which is zero for a product state, and unity for a maxi- 
mally entangled state. This measure is equivalent to the 
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Renyi entropy 12 and is not well suited for much more 
than to single out a pure state, as with increasing n any 
entangled state will converge to zero in this measure. 

In figure 7 we plotted the entanglement according 
to a high order linear entropy and compared with the 
von Neumann entropy. One, possibly not very surpris- 
ing, principal observation is that the numerical calcula- 
tions corresponding to the same parameters, which have 
been considered in figures 4 and 5, gives nearly the same 
behavior (see figure 7). This means that both the en- 
tanglement due to the quantum field entropy Sf and 
high order linear entropy Sl measures are qualitatively 
the same. The important consequence of this observa- 
tion is that the use of quantum entropy and linear en- 
tropy as measures of entanglement are equivalent. We 
must stress, however, that no single measure alone is 
enough to quantify the entanglement in a multi-partite 
systems. Finally, we may say that, it is possible to obtain 
a long living entanglement using the time-independent 
modulating function in pair cat state. This result is quit 
surprising, since the previous studies in the entangle- 
ment for the time-independent interaction for the initial 
coherent state does not contain this interesting feature. 
Which means that the pair cat state as an initial state 
of the field plays very important roles in the quantum 
information process. 

A proposal for the generation of motional pair cat 
states in a two-dimensional anisotropic trap has been 
introduced in Ref. |41| . In order to generate the pair cat 
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Fig. 7 Time evolution of entanglement degree Sl as a func- 
tion of the scaled time At. The parameters are the same as in 
figure 3, where n = 2 (the solid curve) and n = 3 (the dotted 
curve) . 

states for the motion in a two-dimensional anisotropic 
trap we require two laser beams and the frequencies of 
the laser beams are chosen as uoi — ldq — 2vx — '2vy, and 
UJ2 ^ loq. The generation method requires a sequence of 
many measurements of the internal electronic state at 
different time instants appropriately chosen |2j. Apart 
from detecting quantum entanglement in these states, 
might be a useful way to fully characterize simple quan- 
tum logic gates 02 • 



The main results of this work are summarized as follows: 
After using the Heisenberg approach to the quantum 
operator, we have investigated effects produced by the 
ion-field interaction. We have shown that the input pair 
cat states exhibit strong nonclassical features due to the 
correlation between the two modes. The phenomenon 
of quantum revivals in the time-dependent ion-field in- 
teraction and the atomic inversion properties as well as 
the entanglement have been analyzed. Also, it has been 
shown that the atomic inversion as a function of the 
scaled time may display different structures of beats de- 
pending on the initial state of the field and the vibra- 
tional motion as well as on the time-dependent interac- 
tion. Effects such as suppression or attenuation of the 
Rabi oscillations and long time scale revivals as well as 
a periodic dynamics have been observed. 

Using pair cat states as an initial state of the fields, 
our analysis shows that the fluctuations of the entangle- 
ment between maximum and minimum values are irreg- 
ular for both the time-independent and time-dependent 
cases at a short period of the interaction time. Our re- 
sults reveal a new feature in the entanglement as a func- 
tion of the scaled time, that is a long living entangle- 
ment. This feature depends on the parameters ^ and q. 
These interesting features may pave the way to quantum 
information applications such as deterministic all-optical 
quantum computation. The results can be relevant for 
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basic studies on entanglement as well as for applications 
in the preparation and manipulation of quantum states 
for quantum information purposes. As far as scalability 
is concerned, one could think in having several trapped 
ions in a high-Q cavity, exchanging information via the 
cavity photons, or other more sophisticated schemes in- 
volving trapped ions, phonons and photons. 

It is to be remarked that, the entropy of entangle- 
ment can be determined indirectly by experiment us- 
ing the quantum state tomography jllj . This fact comes 
from the definition of the entropy of entanglement which 
is associated with the density matrix of quantum state. 
We believe that the pair cat states may be a useful tool to 
be considered in obtaining a long living entanglement of 
systems of indistinguishable particles. While our analysis 
was carried out with respect to the example of a single 
trapped ion beyond Lamb-Dicke regime, we emphasize 
that the results here reported apply equally well to any 
bipartite system. We hope to report on such issues in a 
forthcoming paper. 
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